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(N . ABSTRACT 



X-ray activity of Anomalous X-ray Pulsars and Soft Gamma-Ray Repeaters may re- 
sult from the heating of their magnetic corona by direct currents dissipated by magnetic 
reconnection. We investigate the possibility that X-ray flares and bursts observed from 
AXPs and SGRs result from magnetospheric reconnection events initiated by devel- 
opment of tearing mode in magnetically-dominated relativistic plasma. We formulate 
equations of resistive force-free electrodynamics, discuss its relation to ideal electrody- 
(T) , namics, and give examples of both ideal and resistive equilibria. Resistive force-free 

current layers are unstable toward the development of small-scale current sheets where 
£f~) \ resistive effects become important. Thin current sheets are found to be unstable due to 
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the development of resistive force-free tearing mode. The growth rate of tearing mode 



is intermediate between the short Alfven time scale ta and a long resistive time scale 



772: r ~ ^/(trta) 1 / 2 , similar to the case of non-relativistic non-force-free plasma. We 



^ ■ propose that growth of tearing mode is related to the typical rise time of flares, ~ 10 

msec. Finally, we discuss how reconnection may explain other magnetar phenomena 
and ways to test the model. 
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1. Introduction 

Two closely related classes of young neutron stars - Anomalous X-ray Pulsars (AXPs) and 
the Soft Gamma-ray Repeaters (SGRs) - both show X-ray flares and, once localized, quiescent 
X-ray emission (Kouveliotou et al. 1998, Gavriil et al. 2002; for recent reviews see Mereghetti 2000, 
Thompson 2001). The energy powering the X-ray luminosity in these sources is supplied by the 
dissipation of super-strong magnetic fields, B > 10 15 G (Thompson & Duncan 1996). Hence the 
two classes are commonly referred to as magnetars. Magnetic fields are formed through a dynamo 
action during supernova collapse (Duncan & Thompson 1992, Thompson & Murray 2001). 

A number of evidence suggest that processes that lead to the production of X-ray flares (and 
possibly of the persistent emission) on magnetars is similar to those operating in Solar corona (see 
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also Section 6; for alternative model addressing persistent emission see Heyl & Hernquist 1998). The 
bursting activity of SGRs is strongly intermittent (good statistics for the AXPs' bursting properties 
does not exists yet). The studies of statistics of SGR bursts from SGR 1900+14 (Gogiis et al. 1999) 
have found a power law dependence of the number of flares on their energy, dN/dE ~ E a , with 
a = 1.66; this is similar to solar flares, where a = 1.5 — 1.7 (Aschwanden et al. 2001). The 
distribution of time intervals between successive bursts from SGR 1900+14 is consistent with a log- 
normal distribution, also similar to the Sun (Gogiis et al. 1999). In addition, Woods et al. (2001) 
have argued that the magnetic field of the neutron star in SGR 1900+14 was significantly altered 
(perhaps globally) during the giant flare. Given these similarities, a suggestion that magnetar 
bursts originate in the current-carrying magnetospheres is only natural. 

Thompson, Lyutikov and Kulkarni (2002) have investigated the global structure of neutron 
star magnetospheres threaded by large-scale electrical currents. In the magnetar model for the 
Soft Gamma Repeaters and Anomalous X-ray Pulsars, these currents are maintained by magnetic 
stresses acting deep inside the star, which may generate both sudden crustal fractures and more 
gradual plastic deformations of the rigid crust. They showed that dissipation of the internal (to the 
neutron star) twisted magnetic field may be done efficiently in much worse conducting magneto- 
sphere where currents are pushed out by electro-magnetic torques. In addition, a significant optical 
depth to resonant cyclotron scattering is generated by the current carriers. Resonant scattering 
in the magnetosphere generates non-thermal component through Compton effect and modifies the 
pulse profile. 

Several possible processes can lead to an explosive release of magnetic energy in magnetosphere 
(Thompson & Duncan 1995, Thompson et al. 2002). First, a sudden twist may be implanted into 
the magnetosphere due to unwinding of the internal magnetic field. This process is accompanied by 
a large-scale displacement of the crust - "gated" by the crust fracture properties. This process is 
likely to occur in a more brittle (as opposed to plastic) crusts. Alternatively, in close analogy with 
Solar flares, a slow, plastic motion of the crust implants a twist (current) in the magnetosphere on 
a long time scale. At some point a global system of magnetospheric currents and sheared magnetic 
fields loses equilibrium and produces a flare. This mechanism has the advantage that the energy 
stored in the external twist need not be limited by the tensile strength of the crust, but instead 
by the total external magnetic field energy. Both these mechanisms advertise that the magnetic 
energy is released in the strongly magnetized plasma through some kind of reconnection. In this 
paper we investigate the underlying electro-magnetic processes responsible for this dissipation. 
The two mechanisms for the production of flares - crust quake and loss of magnetic equilibrium 
- are expected to have many similarities since they both involve dissipation of magnetic energy in 
the magnetosphere. Thompson et al. (2002) and Lyutikov (2002) have discussed how these two 
possibilities may be distinguished observationally (see also Section 6). 

In our opinion the latter mechanism of flare production (loss of magnetic equilibrium) is more 
promising. Thus, we imagine that the magnetosphere of magnetars are similar to the solar corona: 
large multipoles contribute considerably to the total surface magnetic fields; new current-carrying 
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magnetic flux tubes rise almost continuously into the magnetosphere gated by rotational deforma- 
tions of the neutron star crust; as a result, magnetosphere consists of a complicated network of 
interacting current-carrying magnetic flux tubes; coronal fields respond to slow injection of mag- 
netic flux and currents by evolving through a series of quasi-static equilibria, which at some point 
become unstable to resistive reconnection; during relaxation some of the magnetic energy associ- 
ated with the current is dissipated; at a given time magnetosphere may have several active regions, 
where the flux emergence is especially active; dissipation of currents occurs on a wide variety of 
spatial and temporal scales in a very intermittent fashion, with occasional giant flares which involve 
the restructuring of the whole magnetosphere. Overall, heating of the corona and production of 
flares is done by direct currents, which are dissipated by magnetic reconnection (e.g. Browning & 
Priest 1986). 

One of the arguments raised against reconnection is that the typical time scales for the burst 
on-set, ~ 10 msec, is much longer than the Alfven travel time through magnetosphere, which is 
of the order of the light travel time, ~ 30/xsec (Thompson, private communication). In this paper 
we explore a possibility that the flare on-set is due to the development of a tearing mode in the 
resistive magnetically-dominated magnetosphere and show that the typical growth of the tearing 
mode occurs on a time-scale intermediate between the very short Alfven time scale and a long 
resistive time-scale. Tearing mode is one of the principle unstable resistive modes, which plays the 
main role in various TOKAMAK discharges like sawtooth oscillations and major disruptions (e.g. 
Kadomtsev 1975). Tearing mode is also the principle model for the unsteady reconnection in Solar 
flares (e.g. Shivamoggi 1985, Aschwanden 2002) and Earth magnetotail (e.g. Galeev et al. 1978). 

Reconnection in relativistic magnetically-dominated plasma may be qualitatively different from 
the non-relativistic analogue. In force-free plasma currents are allowed to flow only along the 
magnetic field in the plasma rest frame. It is not obvious that any reconnection may occur at all, 
since in conventional models (e.g. , Sweet-Parker), resistive currents flow mostly across magnetic 
field. One example of a qualitative difference of force-free and non-force-free fields is a X-type 
point. One cannot construct a non-trivial (current-carrying) X-type point configuration of a force- 
free field. Another important modification concerns Ohm's law. In the force-free limit Ohm's law 
simplifies considerably since both the inertia of both signs of changes can be neglected and there 
are no Hall terms, which arise due to different masses of charge carries 2 . This bears important 
consequences for the models of reconnection in relativistic magnetically dominated plasmas. In 
non-relativistic plasmas inertial terms in Ohm's law have been invoked to produce what came to 
be known as collisionless reconnection. Typically collisionless reconnection effects are related to 
ion inertial or cyclotron scales. In force-free plasmas (which may be dominated by pairs) no such 
effects arise. 

In addition, the conditions in magnetospheres of neutrons stars, and especially magnetars, 
are different from the better understood conditions in Solar chromosphere in several other ways 



2 Hall term are also absent in pair plasma, e.g. Blackman and Field (1993) 
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that are likely to influence reconnection. First, radiative cyclotron decays times are extremely short 
thus making particle distribution one-dimensional and forcing the currents to flow exclusively along 
the field lines. Secondly, the typical collisional resistivities, which are far too small to explain the 
short dynamical time scales of reconnection of the Sun, are even more suppressed by super-strong 
magnetic fields (e.g. for k^T <C Hlob the suppression is by a factor ~ (^T /Hujb) 2 )- 

In spite of the obvious limitations of the resistive reconnection models, in particular formally 
long time scales, resistive diffusion is known to be able to drive instabilities with time scales much 
shorter than the resistive time scale tr ~ L 2 /r]. This is done through a formation of a very narrow 
current sheet where both the time scales for diffusion may be short and, in addition, resistivity 
may be enhances due to development of plasma turbulence (anomalous resistivity). The resulting 
current sheet tend to be unstable to transverse, k • B = 0, perturbations. This resistive instability 
is called a tearing mode. 

Tearing mode is quite complicated (Furth et al. 1963, White 1983, Fig. 3). The basic phys- 
ical processes leading to development of tearing mode are the following. A current sheet may be 
represented as a set of elementary current filaments. Since current filaments attract each other the 
current sheet is unstable against pairing of current filaments. For sufficiently narrow current sheets 
the energy released by pairing of currents is more than the energy spent to maintain perturbed 
magnetic field - tearing mode has negative energy, so that its dissipation should increase its am- 
plitude leading to dissipation. Magnetic islands form in plasma as a result of the development of 
tearing mode. Typical time scale for the development of the tearing mode is intermediate between 
a short Alfven and long resistive time scale. 

2. Resistive force-free electrodynamics 

Under force- free approximation it is assumed that the plasma dynamics is completely controlled 
by magnetic field. The validity of this approximation requires b 2 3> pc 2 , where b and p are plasma 
rest-frame magnetic field and density. It is assumed that the plasma provides currents and charge 
densities required by the dynamics of electro-magnetic fields, but these current carry no inertia. 

In addition to being magnetically-dominated, microscopic plasma processes, like particle colli- 
sions or plasma turbulence, may contribute to resistivity and thus make plasma non-ideal. Resis- 
tivity will result in the decay of currents supporting the magnetic field; this, in tern, will influence 
the plasma dynamics. We wish to explore how the decay of inertialess currents will affect the 
dynamics of the system. We assume that plasma resistivity can be represented by phenomeno- 
logical parameter r]. To calculate r\ from microscopical principles one needs to take the particle 
dynamics into account and thus go beyond the inertialess approximation. The answer will depend 
on the types of current-carrying particle (e.g. electron-ion or pair plasma) and on the dominant 
type of scattering (e.g. particle-particle, particle- wave or particle- wave-particle) . We wish to avoid 
this complication by introducing a macroscopic parameter r\. We also assume that the dissipated 
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energy of the magnetic field leaves the system, e.g. in a form of radiation. If this was not the case, 
the particle pressure would build up to equipartition breaking the force-free assumption. 

Strictly speaking, inclusion of resistivity violates the force-free condition since non-zero re- 
sistance implies that there is extra force acting on charge carries. Still, in strongly magnetized 
plasmas the resistive forces act only along the direction of the magnetic field, while the dynamics 
transverse to the magnetic field remains unaffected, given by the force-free conditions (12). What 
distinguishes resistive and ideal cases is the way the current is related to the field (Ohm's law). It 
is in this sense that we will use the term resistive force-free electro-dynamics. 

The procedure described above is not entirely consistent in case of relativistic plasmas. The 
reason is that in an ideal force-free plasma the velocity along the field is not defined. Since plasma 
resistivity must be defined in the plasma rest-frame this creates a principal ambiguity. Note, 
that in case of non-relativistic force-free plasma such ambiguity does not appear: it is possible to 
write a self-consistent system of equations describing evolution of non-relativistic resistive force- free 
plasma (Chandrasekhar &; Kendall 1957, Low 1973). In the non-relativistic plasma evolution of 
resistive force- free fields proceeds in a self-similar ways: magnetic fields keep their configuration 
while changing in magnitude. In a process of doing so small electric fields and charge densities 
develop, whose dynamical contribution is neglected in the non-relativistic case. In relativistic case 
we cannot neglect these effects so that evolution force-free resistive fields becomes much more 
complicated. 



Relativistic form of Ohm's law in magnetized plasma may be formally written assuming that 
in the plasma rest frame the current is proportional to the electric field. Several attempts have 
been made to derive Ohm's law in relativistic plasma; none seem completely satisfactory so far. 
Blackman and Field (1993) assumed that there exists a frame where charge density, momentum 
density and velocity flux density are all zero, but in a general case the plasma may be charged 
in its rest frame (defined, for example, as a frame where momentum density is zero). Gedalin 
(1996) wrote expressly covariant form of the relativistic Ohm's law in isotropic plasma. In strong 
magnetic field we expect that transverse and parallel (with respect to magnetic field) conductivities 
are different, so that conductivity becomes a tensor: 



where rf is a conductivity tensor. To find an expression for rf we need to consider particle 
dynamics taking into account particle-particle interaction. Derivation of a general expression for 
rf k is a project worth a separate paper. Here we give simplified expressions for rf k in case of 
strongly magnetized plasma. 



2.1. Ohm's law in magnetically- dominated plasma 




(1) 



To derive Ohm's law in force-free plasma, which assumes a one-fluid description, we need to 
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define the plasma rest-frame. In relativistic plasma the choice of plasma rest-frame is not unique 
(e.g. de Groot et al. 1980). In a force-free plasma the ambiguity of the rest frame is partially 
removed: though the motion along the magnetic field still remains undefined, the motion of particles 
across the magnetic field consists only of an electric drift 

v ^ = — > ( 2 ) 

which allows one to introduce a plasma frame (up to a boost along the field), as a frame where 
drift velocity vanishes. Since both species drift with the same velocity, there is no relative motion 
of particles across the field and thus there is no resistivity across the magnetic field (the current 
across the field is due to exclusively different charge densities of two species). An alternative way 
of reasoning supporting this conclusion is that in force-free plasma there is no current across the 
field in the plasma rest-frame, where resistivity should be defined. 

Thus, under force- free approximation the resistive effects affect only the motion along the field. 
To derive the corresponding Ohm's law we introduce a four-dimensional magnetic field vector 

h = *F ik U k (3) 

where *Fik is a dual electro-magnetic tensor. If B denotes the plasma rest frame magnetic field, 
then Bi = {0, B} and 

6 l = 7 {(v-B),B} (4) 

Convolving eq. (1) with bi we find 

F k U k * F U U 1 = rt% (4 - U k U l ) j t (5) 

Since we expect that in the strongly magnetized plasma resistivity is important only along the field, 
we assume 

rf k = ^ (6) 

We find then 

F lk U k * F u U l = ri(jb) = r, * F a U l f (7) 
since (bU) = 0. Equation (7) gives Ohm's law in force-free plasma. 3 In 3-D notations this gives 

7 2 ((E • v)(B • v) - (E + v x B) • (B + v x E)) = 7/7 (j (B • v) - (J • B)) (8) 

Next we separate electric fields and velocity into components along and across the magnetic field, 
E = E|| + Ej_, v = V|| + v ± with v ± = E ± x B/B 2 . Then eq. (8) simplifies 

E ll (vl-l)=^-(j v ll -J ll )) (9) 



3 Another possibility to define Lorentz-invariant Ohm's law in force-free plasma (suggested by Blandford, private 
communication) is to postulate that r)(jb) = F lk * Fit = E ■ B. The two definitions are different by the term involving 
(B-v)(E-v) in (8). 
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Since the velocity along the field cannot be specified within the framework of force-free plasma 
we chose vu = 0. The ambiguity for the choice of vu stems from the fact that the force-free 
approximation cannot in principle describe the plasma dynamics along the magnetic fields - however 
strong the field is, it does not affect particle motion align the field, so that one must use full MHD 
equations. For example, in a frame-work of two-fluid hydrodynamics there are two equations of 
motion describing the dynamics of each species along the field. The difference of these equations 
gives the generalized Ohm's law. The sum gives the equation of bulk motion of plasma due to the 
influence of external fields. Since we neglect the Poynting flux associated with the bulk motion 
this equation is neglected. Then the choice of parallel velocity becomes just a choice of coordinate 
system. 

In spite of limitations concerning the parallel dynamics of plasma there is a number of relevant 
problems where one can neglect the plasma dynamics along the field either because the variation 
in that direction are small or because of a symmetry of the system (e.g. if magnetic field is directed 
along a cyclic variable). 



Setting v» = in (9) we find 



(E-B) = ^(J-B) (10) 

7 



The factor 7 = 1/ — (Ex B/i? 2 ) 2 comes from the Lorenz transformation of the electric field to 
the plasma rest frame (Jm is invariant under Lorentz transformation along the direction orthogonal 
to magnetic field). Reinstating the transverse current, we find Ohm's law in relativistic force- free 
electro-dynamics 

J _ divj^B) + 1 (E_B)B 

Taking a vector product with B in Ohm's law we find a cross-field dynamical equation 

J x B + divE - ^^ B^j = J x B + E ± divE = (12) 
This generalizes the force-free condition for resistive plasma. 



2.2. Equations of resistive force-free electro-dynamics 

The resistive force-free electrodynamics may be derived from RMHD in the limit of vanishing 
plasma inertia w, p, p — ► and by reinstating electric field 

(E-B)B , . 

E = E|| +E ± = v ; vxB. (13) 

In a force-free formulation the plasma velocity is defined only up to an arbitrary Lorenz boost along 
the magnetic field direction. The velocity across the magnetic field is just the electro-magnetic drift 
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velocity E x B/i? 2 . The dynamical equations are Maxwell equations and Ohm's law (11) 

B = -curlE 

E = curlB -J (14) 

(we use a system of units with the speed of light set to unity; we also incorporate the coefficient 
4-7T into definitions of currents and charge densities). 

Resistive force-free electrodynamics is more complicated than its non-relativistic analogue. In 
the non-relativistic case currents can be expressed from the second Maxwell equation and electric 
field can be eliminated from Ohm's law. In relativistic case this is impossible - the best that we 
can do is to eliminate the current from Ohm's law and solve Maxwell equations for a B and E 
(this is a system of 5 equations since the constraint divB = must be satisfied at all times). 
Resistive force-free electrodynamics is also more complicated than the ideal one. In ideal force-free 
electrodynamics the condition E • B = reduced a number of equations to 4 (Komissarov 2002, 
Lyutikov & Blandford 2003). In addition, in ideal case the second electro-magnetic invariant is 
always larger than 0, B 2 — E 2 > 0, which implies that there is a reference frame where electric field 
is equal to 0. In case of resistive electrodynamics B 2 — E 2 may be smaller than zero, which implies 
an electric field in the plasma rest frame. 

The form of Ohm's law (11) may not be appropriate for some applications: for example, the 
transition to the ideal case (77 — > 0) leave us with no conditions on the current. To make the 
transitions from the resistive to the ideal case clearer we combine Maxwell equations by taking 
scalar products with B and E 

(J • B) = (B • V x B - E • V x E)B - d t (E • B) (15) 

The term B • V x B is called current helicity, the corresponding term involving electric fields does 
not have a name (it does not appear in non-relativistic analysis). 

Using the Ohm's law (11) we can eliminate (E • B) from (15) 



(J • B) + d t [vy 1 - ((E x B)/B 2 ) 2 (J • B)J = (B • V x B - E • V x E)B (16) 

which gives an alternative form of Ohm's law. This form has an advantage that, unlike eq. (11), it 
clearly shows how to make a transition to the ideal case. 

Given Ohm's law (11) we can write down the basic law of energy and momentum conservation 
in resistive force-free electro-dynamics: 

\d t (B 2 + E 2 ) = divE x B — 1 ^ ^ 



d t (E x B) = curlB x B - curlE x E - J x B 

curlB x B - curlE x E + divE (b - j (17) 
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In ideal plasma, rj = 0, Ohm's law gives E • B = 0. This condition plus the Maxwell equations 
(14) are then sufficient to express the current in terms of fields, 

(E x B)V • E + (B • V x B - E • V x E)B 

3 ~ 7^2 ( 18 ) 

while eq. (12) gives 

JxB + EdivE = (19) 

Equations (14) and (18) (or (19)) are then the equations of ideal force-free electrodynamics (Uchida 
1997, Gruzinov 1999, Komissarov 2002, Lyutikov & Blandford 2003). 



2.3. Applicability of force-free approximation 

Force-free electro-dynamics assumes that inertia of plasma is negligible. This approximation 
is bound to break down for very large effective plasma velocities, when electric field becomes too 
close in value to magnetic field, E — > B. The condition that inertia is negligible is equivalent to 
the condition that the effective plasma four- velocity, u ~ E/(By/l — (E/B) 2 ), is smaller that the 
Alfven four-velocity in plasma. Assume that in the plasma rest frame the ratio of magnetic energy 
density to plasma energy density (which may include both rest mass energy density and internal 
energy) is a. In a strongly magnetized plasma d > 1 (in a force-free plasma o is assumed to be 
infinite). Then the Alfven wave phase velocity and the corresponding Lorentz factor are (e.g. 
Kennel & Coroniti 1984, Lyutikov & Blandford 2003) 

u\ = a (20) 

This puts an upper limit on the value of electric fields consistent with force-free approximation 

B 2 -E 2 



E 2 »* < 21 > 



3. Relativistic Current Sheets 

3.1. Ideal Force-free Current Sheets 

Since we are interested in the stability of weakly resistive configurations, we first consider 
steady state planar solution of the ideal force- free electrodynamics (eqns. (14) and (18)). Assume 
that magnetic field lays in the x — y plane, B z = 0, and that all quantities dependent only on z. 
Then from the first Maxwell equation we find that E x , E y ~ const. The x and z components of 
the curl B = J equation then give (y component is the same as x) 



d z {Bl + Bl-El)=V 
(B y E x — B x E y ) d z E z = 



(22) 
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In addition, for ideal plasma condition E • B = requires 

B X E X + B y Ey = (23) 

Below we discuss several possible solutions of these equations. The two main types of current 
sheet equilibrium are the current sheets (when the electric field is zero and there is no motion of 
plasma) and shear layer, where non-zero electric fields lead to plasma motion in the plane of the 
current sheet or into it, creating magnetic and velocity-sheared configurations. 

1. Crossed constant electric and magnetic fields. We assume that the components of the electric 
field in the plane of the sheets are non-zero, E x E y / and, in addition, E z ~ const. This 
case corresponds to plasma in constant electric and magnetic fields moving with velocities 



V 



E z E y E z Eg Ey 



B^l + El/El B E x Jl + EH El ' B E xy /l + E%/E 



(24) 



(in Cartesian coordinates x, y, z). By Lorentz transformation electric fields can be completely 
eliminated, so that in the plasma rest frame E = 0. 

2. Sheared force-free layer. Assuming E x = E y = 0, E z / we find 

Bl + Bl-E 2 z =Bl (25) 

where Bq is some constant. There are two interesting types of equilibrium here. First is a 
sheared B x field in a constant guiding field B y = Bq/Uq (Uq is some constant). 

B = B {tanh z/L, -^,0 
E = 5 {0,0, tanh z/L} 
J = -£{0, sech 2 z/L,0} 

p e = — y~ sechz/Ltanhz/L 

Y-{ t/otanhz/L Uq tanh 2 z/L 

~ 1 1 + C/ 2 tanh 2 z/L'~T+ C/ 2 tanh 2 z/L ' 

2 C/n tanh 2 z/L 2 ^ a . 2 lr 

v 2 = 2-= 4—-, 7 2 = 1 + Uq tanh 2 z L 

1 + C/ 2 tanh 2 z/L ' ' 

P = B 2 |^tanhz/L,tanh 2 z/L| (26) 

where L is a characteristic width of the layer, P is a Poynting flux. Thus, there is a Poynting 
flux of electro-magnetic energy along the current layer. At the different side of the current 
layer P y changes sing, while P x remains the same. (Fig. 1). The layer is sheared both in 
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x and y direction. On the different sides of the current layer the plasma is streaming in 
opposite direction along the x, with a four velocity at infinity reaching u x = ±Uo/y/l + Uq , 
while along the y direction the plasma is streaming in the same direction on both sides of the 



layer, reaching at infinity u y = —Uq/^JI + Uq. 

Another possible type of layer is a sheared and rotating magnetic field, (Blandford, private 
communication) 



B 



Bq |tanh z/L, \[2 sech z/L, o| 



E = S {0,0, sech z/L} 

B ( i 

— |\/2 sech z/L tanh z/L, sech 2 z/L,0 j 



J 



p e = — p sech z/L tanh z/L 



v = J sinh z/L Q 

' 1 + cosh 2 z/L ' 1 + cosh 2 z/L ' 



v 2 = —s — — , 7 2 = 1 + sech 2 z/L 

1 + cosh 2 z/L ' 



|-\/2sech 2 z/L, sech z/L tanh z/L, o} (27) 



(Fig. 2). 



3. Magnetic rotation discontinuity E = 0. In this case we find B 2 = B 2 , so that magnetic field 
rotates keeping its absolute value constant (Fig. 3). This case will be our primary interest. 

Finally we note, that since ideal electro-dynamics supports only two types of waves (luminal 
fast modes and subluminal Alfven modes), the current sheets discussed above are nothing else 
but Alfven waves considered in some preferred frame of reference (Komissarov 2002). Then the 
instability considered below are related to the instabilities of 1-D dissipative Alfven waves (we 
thank S. Komissarov for pointing this out). 



3.2. Resistive Force-free Current Sheets 

In a strongly resistive media one may seek stationary solutions of the equation of resistive 
electro-dynamics (14, 11), assuming that resistive terms are important everywhere (in a way similar 
to important works of Low (1973) for non-relativistic plasma). Assuming a sheared magnetic field 
field configuration with E z = we find 

B > = 1 B X (E • B) 

y V y/B* - (E y B x - B y E x y 

1 BJB-B) , N 

K = ; g 28 

V y/B* - (E y B x - B y E x y 
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From which it follows that B 2 = const. Assuming B x = cos <j){z), B y = sin^(z) and setting without 
loss of generality E y = we find equation for cf>: 

dA = ~- -£^£= (29) 
V v 7 ! -El sin 2 4> V ' 

Which can be integrated to give an implicit dependent 4>{z): 



E x z E x . E x . y/l - (E x j Bq) 2 

= — — arcsm(— simp) — - In .4 

B B 2 




E x \ 2 . , I {Ei 



sin-o • 2^ 1 sin 0yi_^_£j sin^ | |;!()) 

It describes a inflow of plasma into resistive current layer. For small (E x /Bq) — > this gives 

4> = - tanh (31) 
BoV 

which corresponds to the Low (1973) solution 

„ r, , E x z 
B x = B sech — — 

-Bo V 

E x z 

B v = — Bq tanh — — 
BoV 

v * = -jr tanh ir^ (32) 

■DO -DO T/ 

Application of this solution to astrophysical plasmas suffers from the same problem as the 
original Low's solution: the diffusion time scales are much longer than the time scales of interest. We 
suspect (though did not show it) that the relativistic solutions (30) is unstable to small perturbations 
which would lead to electric current singularity. This brings us to the main part of the paper - 
development of resistive current sheets in a relativistic force-free plasma. 



4. Tearing mode stability of rotational discontinuities 

A standard method in describing the evolution of tearing mode is similar to the boundary layer 
problem 4 . It involves a separation of a current layer into a "bulk", where derivatives are small and 
resistivity is not important, and a narrow "boundary layer" where derivatives and resistivity may 
be large. Two different approximations are done in each layer - ideal and weakly varying plasma in 
the bulk and a narrow resistive sublayer (Fig. 3). Two solutions should be matched continuously. 



4 Kinetic approach to tearing mode is described in in Galeev (1984). It was generalized to relativistic particle 
dynamics by Zeleny & Krasnoslskih (1979). 
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In this section we investigate a resistive stability of a relativistic current layer. We assume that 
initially the electric field is zero, so that the magnetic field while remaining constant in magnitude 
rotates over some angle (Section 3.1), which we chose to be ir radians for simplicity. 

Assume that a current layer has a width L. Since B 2 = const, a possible form of the unper- 
turbed magnetic field is 




Next we investigate stability of such current layer to small perturbations. 



4.1. Stability of resistive force- free current sheet 

In this main section of the paper we investigate the stability properties of a planar resistive 
current sheet. We assume that in the bulks of the plasma resistivity is small, so that initial 
configuration is described by an ideal current sheet in which magnetic field remains planar {B z = 0) 
and rotates over some angle B 2 = B 2 , while electric field is zero E = 0. Formally we should start 
with resistive solution (30) with has non-zero velocity of plasma into the current layer, but for large 
scale current that resistive inflow velocity is negligible. 

Consider small fluctuations of the electro-magnetic fields and current 

B = B + b 
E = e 

J = Jo+j (34) 

Assume that the perturbations vary as exp{i(wi — k x x — k y y)}. Then eqns (14), (16) and (12) give 

b z uj — (k x e) z = 

i(k-b) -d z (k x e) 2 = 

e z k 2 + uj(k x h) z u — id z (k ■ e) = 

(k-j) + iu(k-e) + d z (k xb) 2 = 

b z k 2 + i(k x j) 2 - w(k x e) 2 - id z (k • b) = 

j z + i({k x b) z + ue z ) = 

(B o -e)-»/((Bo-j)-(Bo-b) J ff) = 

j z + b z H = (35) 

where H = (B x d z B y — B y d z B x )/B 2 is a normalized current vorticity of the initial state and 
k = {k x , ky , 0}. 
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We can readily eliminate current, (k x j) z , (k • j), j z , as well as z-components of electric and 
magnetic field fluctuations e z and b z . After some rearrangements the equations involving fluctua- 
tions of x- and y— components of electric and magnetic fields become 

i(k-b)-d 2 (kxe) 2 = 

u 2 \ „ „ , to 2 



ico(k x b) 2 ( 1 - -p-J - Hu(k x e) 2 - ^^(k • e) = 

(B ■e)-r t (i(k x B ) 2 (k x e) 2 - 2luH(B • b) - ud z (B x b) 2 ) = 

i(k • B )(k x e) 2 - iuj 2 (B x e) 2 + d 2 (B • b) = (36) 

From which we can eliminate magnetic fields 

m i_\ ■ /(k-B )a 2 (kxe) 2 (kxB ) 2 (K 2 (kxe) 2 +A(k.e))\ 

( B 0- b ) = ^ ^ ^ + fcMfc2- W 2) ) 

m i •> ./ (kxBo),3,(kxe), (k • B ) (g (k x e) 2 + u; 2 /fc 2 g 2 (k • e)) \ 

(B ° X b)z = ? ^ fc^ u,(fc 2 -c 2 ) ) (37) 

This, with corresponding equations for electric fields 

(B • e) - n (t(k x B ) 2 (k x e) 2 - 2#cu(B • b) - iu 2 (B • e) - ojd z (B x b) 2 ) = 

i(k • B )(k x e) 2 - iw 2 (B x e) 2 + d z (B • b) = (38) 

is where we stop the general approach. 



4.2. Ideal current layer 

In the bulk of the plasma we can neglect resistivity so that the dynamics of the current layer is 
governed by equations (14,18). In this section we investigate stability of an ideal force-free current 
sheet to small perturbation (see also Gruzinov 1999). 

For ideal plasma i] = 0, we find 

d z (Fd z -^-) + -^F(u 2 -k 2 ) = (39) 
Or, returning to the definition of b z 

where F = (B • k) 2 — uj 2 Bq. Equation (40) describes stability properties of an ideal relativistic 
force-free current sheet. A quick examination shows that similar to the non-relativistic case points 
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(Bo • k) = is a special point of the equation. Fluctuation of the magnetic field b z diverge at these 
surfaces creating current sheets. 

Consider fluctuations near the surfaces (B • k) = 0. Assume that B x becomes zero at z = 0; 
then, neglecting variations in y direction, k y — > and introducing displacement £ = b z /B x k x we 
find 

[Blkl - Blu?) {$(u? - kl) + d 2 z i) + 2B x tf x d£d z B x (41) 
which can be further simplified 

{kl - Cj 2 ) {Blkl - Cj 2 ) t(Z) - dim (42) 

where 



k x — kx L 
j = 

B x 



Cj = uL 
B x 

Bo 



z 

Z= L 



I 



k 2 x B 2 x - Cj 2 



dz (43) 



For example, for a linear B x = z, we find Z = —avctan(k x z /Cj) / (k x Cj) . Then eq. (42) gives 

Cj\k 2 x - Cj 2 ) sech 2 (k x CjZ)t - 8U{Z) = (44) 

which has a form of a non-linear oscillator. 

General solution of Eq. (44) is quite complicated. Its dependence on Cj shows that variations 
of the field have electro-magnetic structure (as oppose to purely magnetostatic variations in the 
non-relativistic case). We expect that typical growth rates will be small Cj <C k x , so that we can 
neglect the electro-magnetic corrections by setting Cj — > 0. This will have to be checked o posteriori. 
This limit corresponds to non-relativistic approximation, so that the relations (45-48) coincide with 
the familiar relations for the non-relativistic tearing mode. 

In the limit Cj — > eq. (41) gives for B x = z when \z\ < 1 and B x = 1 for \z\ > 1 

[ (~i coshfctti , n. smh k x z -f | - i ^ i 
£=) Gl z 1 II \Z\ < 1 ^ 

i Co exp{— if \z\ > 1 

where we have chosen a solution which is decaying as \z\ — > oo. Matching the two solutions at 
\z\ = 1 gives 



d = 5°- (2fc x - 1 + exp{-2fc B }) 
Co 



C 2 



f 1 - 2^ + exp{-2fc B }) (46) 



2k x 
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Given a solution for £ we can find magnetic field fluctuation 

h 



- = B x k x £, = k x I C\ cosh k x z + C2 sinh k x z ) (47) 
b ^ ' 

Magnetic field is continuous at z = but its derivative experiences a jump 

2C2k x 2 



A = 



b 



z 



— , for ^ < 1 (48) 



This means that there is current sheet forming at z = 0. 



4.3. Structure of resistive sublayer 



To study the structure of the resistive sublayer we make several approximations for equations 
(37-38). First, we neglect variations along the y-axis, k y — > 0. We find then 

e x B x + e y B y + ^ (e y (S^fT 2 + k 2 x ) + B x d z H) + 2B x Hd z e y - B y d 2 z e y + 

-jj-^f2 (( B * e * + B y e y)( k l + r ') + + B x 8 z H)e y + H(B y d z e x + £^e y ) + B^e*)) = 

e y {B x (k 2 x - H 2 ) - B y d z H) - B x d 2 e y + 

p2 

k 2 + T 2 ^ e v B * ~ B vex)(*Z + T 2 ) + e y (B 2 H 2 + B y d z H) + H(B x d z e x + Sj,^) + B y d 2 e x ) =(89) 
where T = —iio. 



We expect that the width of the resistive layer, ~ eL is much smaller than the width of the 
current layer L, e <C 1. We then expand system (49) near z ~ e by introducing ( = ez/L. 



For e <C 1 we find 



<9 2 e y 



c(#-^) + ^(&-S<aa))=o (5o> 



hi v « c 

where 77 = t]/(TL 2 ), T = TL, k x = k x L (T = TL/c in dimensional units) and we used B x ~ 
i?o-z, -By ~ -Bo, <9z-B ~ —z/L 2 , H ~ —1/L. We also assumed 77, T <C 1 and neglected terms of higher 
orders in parameters k x ,T],T 2 . 

Eliminating e x we find a fourth order equation for e y (in which we can also neglect (eC 1 and 
f < k x ): 

C f 2 t) (C<9 f % - 2dfe y ) + (2f 2 t) - t 2 e\t 2 + ( 2 k 2 x e 2 )) d 2 e y + 4 C f 2 e 2 d ( e y - e 2 ( 2 f 2 - fc^ 4 ) = 

(51) 
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The typical scales appearing in eq. (51) are 

x ^ 
di = 



e V eTL 2 



\fk~ x t V ek x \fkxL 

5\ is a resistive skin depth (the distance that the field diffuses in time r -1 ), 82 is a width of tearing 
sub-layer, £3 <C 5\ does not seem to have a physical meaning. 

For very small £ — > eq (51) becomes 

<9 f % - = (53) 

which has a solution 



~ cosh — (54) 



\pr\ 

where we chose e' y {C, = 0) = 0. This solution is valid for Q <C 82- For the order-of-magnitude 
estimates we can assume that this solution is valid until £ ~ 82- Setting 82 = 1 we find the width 
of the resistive layer: 

S/L = C = M^ = M (55) 
V ™x V 

(compare with Woods (1987), eq. 7.116). 

Equating then A = e' y /e y taken at ( = 1 to external solution we find 

e 2 e 

A = — = ^ (56) 
V k x 



From which we find 



e ~ f (57) 

fox 



The growth rate then follows from eqns (55) and (57) 

77 c 



3.2 \ 1/5 j 



fc2L 10 J k 2 x L 2 3/5 2/5 



(58) 



The growth rate increases with l/k x . The maximum rate may be found from the condition that 
the resistive time scale for a sublayer of width eL, is much smaller that the growth rate: 

r<(^) 2/3 (^) 1/3 (59) 
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Which gives 



k x L 



\Tc) 



(60) 



and 




(61) 



This estimate formally coincides with the case of non-relativistic non-force-free plasma. 



5. Application to magnetars 



One of main shortcomings of the current approach is that resistivity r] was not calculated from 
the particle kinetics, but was introduced as a macroscopic property of plasma - a common approach 
in continuous mechanics. Resistivity in tenuous astrophysical plasmas is due to collective processes 
and not binary collision. It has to be excited by plasma currents and thus is likely to be a (non- 
linear) function of a current itself. This requires a kinetic treatment of plasma-wave interaction 
as well as a correct account of non-linear (or quasi-linear) feedback of plasma turbulence on the 
particle. Excitation of plasma turbulence by currents in force-free fields is confirmed by a number 
of electro-magnetic particle-in-cell simulations {e.g. Sakai et al. 2001). Qualitatively, when the drift 
velocity of a current exceeds the plasma thermal velocity, strong plasma turbulence is excited. 

In spite of complicated microphysics which determines the resistivity we can make qualitative 
upper estimates on the value of resistivity r]. At the early stages of the development of instability 
the plasma remains force-free, until a large amount of magnetic energy has been dissipated. Since 
under force-free assumption the particles are bound to move only along the field lines, this limits 
considerably a number of possible resonant wave-particle interactions that can lead to development 
of plasma turbulence. The two remaining options are the Langmuir turbulence, which in relativistic 
plasma develops on a typical scale of electron skip depth, 5 e ~ c/uj Pje , and ion sound turbulence, 
which in relativistic plasma develops on an ion skip depth 5i ~ c/uj p ^ (w p , e and to Pj i are the electron 
and ion plasma plasma frequencies) . They are different by a square root of the ratio of electron to 
ion masses 5 e /5i ~ {m e jmi) x l 2 . Qualitatively, a fully developed turbulence with a typical velocity 
c and typical scale 5 would produce a resistivity 



7] 



(62) 



Using this estimate in (61) we find 




(63) 



Note, that the growth rate of the tearing mode is proportional to (m e /m i ) 1 / 4 ', which is only a factor 
of few. 



In current-carrying magnetospheres of magnetars the plasma frequency is (this estimate comes 
from he requirement that toroidal magnetic field created by currents does not exceed poloidal, 
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Thompson et al. 2002) 



^/uj B r/c ~ 2 x 10 9 rad/sec 



(64) 



for B ~ 5 x 10 14 G and r ~ 10 7 cm. 



Then 



5 ~ 10 cm, 7] ~ 3 x 10 11 cm 2 /sec 



(65) 



Resistive and Alfven time scales are then 



t r ~ L 2 /r] = 3Lj 6 sec 

ta ~ L/c = 3 x 10~ 5 L 10 6 sec 



(66) 



where L 10 e = L/10 6 cm. 



The tearing mode growth time is then 



1 




(67) 




For a given L this is the lower estimate on Tt since we used an upper estimate on rj. 

The growth time (67) is of the order of the observed rise time of the SGR X-ray flares, < 10 
msec. For a given current sheet width L eq. (67) gives a lower estimate and also in Since in the 
observed bursts the rise time is limited by the intensity of the burst - weaker bursts are expected 
to have shorter rise times (Gogiis 2002) - smaller burst should come from smaller current sheets. 



In this paper we have first formulated equations of resistive force-free electro-dynamics, which 
describes strongly magnetized dissipative relativistic plasmas. This is an important step towards 
understanding of plasma dynamics in extreme conditions (when energy density is dominated by 
magnetic field energy density) and may serve as a guiding rule for numerical investigations of a wide 
variety of electro-magnetically dominated astrophysical plasmas (since most numerical schemes are 
necessarily dissipative) . 

Next we have analyzed development of unsteady reconnection in a relativistic force-free plasma. 
We find that under assumptions of a force-free resistive plasma the development of a tearing mode 
proceeds qualitatively in a way similar to the non-relativistic case, when inertia of matter is not 
important. This is a bit surprising result, given that we solve a quite different system of equations. 
Several factors may explain this. First, the special relativistic corrections ( — (Ex B) 2 /i? 2 

term in (11) did not enter our linearized analysis since we were considering pure magnetic equi- 
librium, without an equilibrium drift of particles. Thus, the solution for the bulk of the current 
layer virtually coincides with the non-relativistic case, since even in the non-relativistic case the 
motion in the bulk of the current layer is assumed to be inertia- and resistance-free, thus obeying 



6. 



Discussion 
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ideal force-free equations. Secondly, in relativistic regime the important quantity that controls the 
dynamical response of a plasma to a given force is not inertia but energy density. In force-free 
plasmas the energy density is completely composed of magnetic field (one may say that magnetic 
field gives an effective inertia to plasma). 

We find that the typical rise time of SGR flares, ~ 10 msec may be related to the development of 
unsteady reconnection via tearing mode in the strongly magnetized relativistic plasma of magnetar 
magnetospheres. To obtain this estimate we have assumed that the plasma resistivity is provided 
by relativistic Langmuir turbulence with a typical spatial scale of the order of the skin depth, 
~ c/ujp, that the corresponding plasma density is found using the Thompson et al. (2002) model of 
current-carrying magnetosphere, and that the resistivity is given by the fully-developed relativistic 
Langmuir turbulence, rj ~ c 2 uj p . 

The relation between the rise time of a flare and the total energy output cannot be simply 
predicted in this model. On one hand, growth rate is independent of the magnetic field and 
the value of the current in the layer. On the other hand, we expect that the resistivity rj is 
anomalous resistivity excited by currents and thus may correlate with current strength and the 
amount of energy stored in non-potential magnetic fields. In addition, very long burst with multiple 
components should result from numerous avalanche-type reconnection events, as reconnection at 
one point may triggers reconnection at other points (e.g. as in the SOC model of Li & Hamilton 
1991). 

Development of tearing mode is likely to be accompanied by acceleration of particles and 
production of high energy emission similar to the well studied acceleration in the Earth magnetotail 
(e.g. Coroniti and Kennel 1979). Qualitatively, acceleration may be done by electric fields directed 
along magnetic fields, or near magnetic null surfaces. At later stages of reconnection various 
mechanisms of particle acceleration (DC electric fields, stochastic acceleration, shock acceleration) 
may be operational (e.g. Aschwanden 2002). 

An alternative possibility for production of flares on magnetars is that flares may result from 
a sudden change (unwinding) in the internal magnetic field. In this case, a twist is implanted into 
the magnetosphere. A large-scale displacement of the crust probably requires the formation of a 
propagating fracture, close to which the magnetic field is strongly sheared (Thompson & Duncan 
1995, 2001; Woods et al. 2001a). For this scenario it is important that the critical crustal shear 
stress be relatively large, ~ 0.1. It is hard to estimate the critical stress for the neutron star crusts 
(critical stress cannot be calculated from first principles). Indirect evidence based on possible 
measurement of free precession (Cutler et al. 2002) show that critical stress is small, ~ 10~ 5 , 
favoring the plastic creep possibility. 

Below we summarize the evidence that point to the magnetospheric origin of flares (see also 
Thompson et al. 2002). 

• The pulse profile of SGR 1900+14 changed dramatically following the August 27 giant flare, 
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simplifying to a single sinusoidal pulse from 4-5 sub-pulses (Woods et al. 2001a). In the 
reconnection model the post-flare magnetosphere is expected to have a simpler structure, as 
the pre-flare network of currents has been largely dissipated. 

• The persistent spectrum of SGR 1900+14 softened measurably after the 27 August giant 
flare: the best-fit spectral index (pure power law) softened from -1.89 ± 0.06 to -2.20 ± 0.05 
(Woods et al. 1999a). Since the spectral index is a measure of the current strength in the 
magnetosphere (Thompson et al. 2002), this points to weaker currents in the post-flare state, 
consistent with dissipation of currents during a flare. 

• The typical rise-time of bursts, ~ 0.01 sec is consistent with the time-scale for the development 
of force-free tearing mode in the magnetosphere. 

• SGR bursts come at random phases in the pulse profile Palmer (2000) - this is naturally 
explained if (even only one!) emission cite is located high in the magnetosphere, so that we 
see all the bursts (if the bursts were associated with a particular active region on the surface 
of the neutron star, one would expect a correlation with a phase); 

• pulsed fraction increases in the tails of the strong bursts, keeping the pulse profile similar 
to the persistent emission (Woods et al. 2002) - this is easier explained if the energy release 
processes occurring high in the magnetosphere after the giant burst are connected to the same 
hot spot on the surface of the neutron star as the field which are active during the quiescent 
phase. 

• smaller fluency SGR events, have harder spectra than the more intense ones (Gdgiis et al. 
2002) (this is also true for the spikes of multi-structured bursts); this is consistent with short 
events being due to reconnection, while longer events have a large contribution from the 
surface, heated by the precipitating particles. 

We will also make several qualitative suggestions regarding the possible relation between AXPs 
and SGRs in a frame-work of magnetospheric model for energy dissipation. The prime question 
is why AXPs and SGRs look so different (AXPs emit mostly persistent emission, while SGRs are 
actively bursting) in spite of similar physical conditions. First, following the Parker's paradigm for 
Solar flares, it is possible that the persistent emission of magnetars, as well as bursts, is powered 
by small scale reconnection events. The quasi-stable profiles are then due to multiple scattering 
of radiation in the magnetosphere at large radii r ~ 10rs (Thompson et al. 2002). At these radii 
the field is dominated by lower multipoles, which dissipate on longer time scales than the higher 
multipoles. Secondly, AXPs have stronger than SGRs measured magnetic fields and thus can 
support larger non-potential magnetic fields and larger currents flowing through magnetosphere. 
Larger currents create higher level of turbulence which contributes to larger resistivity. In a medium 
with larger resistivity the intermitency is less pronounced, so that dissipation of currents in AXPs 
is dominated by small scale reconnection events, while in SGRs bursts on average contribute as 
much energy as persistent emission. 



- 22 - 



There is a number of relevant problems to be solved. Primarily, there is a need for sim- 
ulations to confirm the growth rates and to test the non-linear stages of the instability. This 
requires development of electro-magnetic codes. We are aware of the ongoing development of ideal 
electro-magnetic codes, but none have been completed so far (MacFadyen, private communication; 
Spitkovsky, private communication). Several exact equilibrium solutions presented in this work 
offer good tests for numerical schemes. Secondly, resistive instabilities of relativistic cylindrical 
plasmas need to be studied as well in connection with the possible role of reconnection in AGN and 
pulsar jets. One expects that in a cylindrical geometry the growth of resistive instabilities will be 
higher than in the planar case. Non-linear development of resistive instabilities is also of interest, 
but given a considerable complications involved in calculating non-linear stages of instabilities and 
obvious difficulties in relating the results to observations, it is not considered a promising approach 
(at least analytically) at the moment. 

I would like to thank Eric Blackman, Roger Blandford, Sergei Komissarov, Anatoly Spitkovsky, 
Vladimir Pariev, Dmitry Uzdensky for their interest in this work. 
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Fig. 1. — First type of electro-magnetic sheared layer - sheared B x field in a constant guiding field. 
At large \z\ — > oo the medium is moving in opposite direction with total four-velocity ±Uq; On the 
different sides of the current layer the plasma is streaming in opposite direction along the x, with 
a four velocity at infinity reaching u x = ±Uo/y/l + Uq, while along the y direction the plasma is 
streaming in the same direction on both sides of the layer, reaching at infinity u y = —Uq/\/1 + Uq. 
The curves, plotted for Uq = 1, correspond to (1) total Lorentz factor 7 , (2) total magnetic field 
B 2 = B x + By, (3) current j y , (4) velocity along y direction v y , (5) velocity along x direction v x , 
(6) B x . 
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Fig. 2. — Second type of electro-magnetic sheared layer - rotating and sheared magnetic field. The 
curves correspond to (1) total Lorentz factor 7 , (2) B y , (3) current j y , (4) velocity along x direction 
v x , (5) velocity along y direction v y , (6) B x . 
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Fig. 3. — Development of a tearing mode. Initial configuration consists of a current layer in which 
magnetic field rotates by it radians. Current is flowing along magnetic field. A thin current sheet 
forms due to development of tearing mode. Magnetic islands form inside the resistive current sheet. 



